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t. INTRODUCTION 

*. I Soma defects of turbulence models 

1» has become customary among engineers and others concerned 
wills practical (low and neat-transfer predictions to use 'turbulence 
models', to sets of equations which purport, when coupled with 
the equations ol mean motion, to describe such statistically 
significant properties as the ‘effective viscosity', 'length scale' and 
'energy*. Most frequently used are those models employing two 
partial differential equations, the dependent variables of which are. 
for example, the turbulence energy, k, and the dissipation rate «£. 

Although valuable service has been performed by such models, 
their predictive ability is far from meeting ail demands. A brief ust 
of their shortcomings Is: 

(a) The 'Iniermtuancy' of turbulent flow, evident to all who took 
perceptively at. for example, automobile-exhaust plumes, is 
left out of account. 

(b) Attempts to extend the models to chemically-reacting flows 
have had Uttle success, because chemical reactions take 
place In regions of steep gradients which find no expression 
in the equations employed. 

<c) Empirical and non-general corrections have had to be made 
to the 'constants* of the equations In order to bring their 
predictions Into tine with experimental evidence regarding the 
influences of mean-streamline curvature and of Interactions 
between gravitational forces and temperature gradients. 

(d> The models are totally unable to explain the 'unmlxlng* 
phenomenon which steepens gradients of average fluid 
properties rather than diminishing them. 

In the view of the present writer, the above defects have a 
common origin, namely neglect of the ‘spotilness* of real turbulent 
flows. The question therefore to be discussed Is: How can the 

major effects of the ‘spottiness* bo described by a mathematical 
apparatus that Is still not too elaborate for practical use? 

1.2 A proposed shift of emphasis 

<a> The origins 

Modern Ideas about turbulence have been greatly influenced by the 
writings of Osborne Beynolds (eg Ref 1) and Ludwig Prandll (eg 
Bel 2> . whose observations had led them to conceive ol turbulence 
as Involving a semi-random exchange of matter, simitar to that 
envisaged In the kinetic theory of gases, but on a larger scale. 
The following quotations Illustrate UUs. 
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Reynoids: 

“The h©at carried off . . is proportional to the rate at which 
particles or molecules pass backwards and forwards from the 
surface*. (Present author's Italics. ) 

PrandU: 

“U Is now necessary to make a usable hypothesis for the 
mixing velocity, w. The transverse momentum associated 
with this velocity must be constantly destroyed by braking. 
and constantly re-established* . (Present author's Italics) 

It Is a consequence of this pre-occupation with mixing that 
turbulent fluxes of a scalar quantity are almost Invariably 
expressed, by analogy with Pick's law of diffusion and Fourier's law 
of heat conduction, in terms of the gradient of the corresponding 
intensive property. Thus, for heat transfer, the taw: 


<q > d " “^eff* 9 r a<J T .(1.2-1) 

is employed, along with the supposition that is a 

greater-than-zero measure of the local turbulent motion. Here 
<q >0 Is the heat-flux vector. x 9 ff is the effective conductivity, and 
T Is the time-mean temperature. 

So embedded is this notion In current thinking that, when 
experimental evidence becomes incontrovertible that the so-defined 
*eff !s negative, the phrase “contra-gradient diffusion* (Ref 3) is 
applied to the- phenomenon, which is regarded as anomalous. 

(to) The sifting phenomenon and Its Implications 

in the view ol the present writer, normal, le “co-gradient* diffusion 
ts only one c» the two main ways in which turouleni fluxes are 
caused: the second, although It could be called “ non -gradient 

diffusion*. Is better described as ’sifting*. This is the 

phenomenon according to which Intermingled fragments of fluid 
differing In. say. density, move relative to one another as a 
consequence of a body-force field. 

The sifting process is most easily seen when surface tension keeps 
the two fluids apart, as when steam bubbles rise from the healed 
bottom of a kettle. This is not an exchange process: for the 

bubbles rise continuously: yet It Is undoubtedly associated with a 

flux of energy. 
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The mathematical expression for the heel flux associated with 
silling is conveniently written as: 


<q> s * constant <a> h' . <1 .2-2) 

which slgntlles that the flux <Q> S depends upon the local 
fluctuations of enthalpy, h*. and on the acceleration vector <a>. 
Similarly, if w-dlrectlon momentum transfer Is In question, the 
contribution ot sifting to the shear-stress vector <t> s may be 
expressed as: 


<r> s = constant <a> w* .(1.2-3) 

a term which may or may not have the same sign as the more 
conventionally defined viscous contribution: 


<T>d “ i*ef.f grad • .u .2-4) 

Here w and w* represent the time-average and fluctuating 
components ot longitudinal velocity. 

tcf Circumstances In which -turbulent sifting" Is significant 
Cl) Bubble rise in a saucepan 

In the phenomenon fust mentioned. It 1$ Important to emphasise, 
the vertical-direction temperature gradient Is very small; for both 
steam and water are at the saturation temperature corresponding to 
the local pressure, which varies negligibly over the smalt height 
change 

Nor can enthalpy or concentration gradients account for the bubble 
motion: for the bubbles are moving towards the region of highest 

steam concentration, namely the top part of the saucepan, not 
away from It. The fluid enthalpy there is ol course higher than at 
the bottom. 

til) Rise of smoke In a room 

if a fire starts on the floor in the corner of a room, flames and 
smoke rise in a turbulent manner, often forming a plume of smoke 
within which the naked eye can distinguish fragments of smoky gas 
separated by less opaque bodies of clear air; yet the layer of 
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smoke which flows along the celling Is usually separated from the 
air below It by a sharp Interface: the lighter smoke fragments 

have floated upwards leaving the denser air behind. 

(til) Boundary layers on curved surfaces 

When a turbulent boundary layer encounters and flows over a 
convex surface, the pressure gradient along the normal to the 
surface is positive. As a consequence, fragments of fluid which 
are flowing at a greater-'han -average velocity are flung outward, 
white slower fragments are forced towards the wall. There Is 
therefore, from this cause. , a momentum transfer outwards, 
towards the region of higher velocity. 

The momentum flow Is reversed when the curved surface Is 
concave: for the higher-velocity fragments move towards the wall. 

The momentum flux associated with sifting has. in this case, the 
same sign as that associated with turbulent mixing, 

Uv> Flow In a cyclone 

Similar sifting flows occur In cyclone separators. That particles of 
denser material move preferentially outwards is of course obvious: 
it Is the raison d'etre of the apparatus. Less well-recognised is 
the fact that fluid particles having larger circumferential velocities 
will migrate tc the larger radii, white the slower-moving particles 
move towards the cyclone axis. 

This steady transfer of kinetic energy to the outside Is presumably 
the reason for the observation that. In the Fanque-HHsch tube 
(which Is a kind of cyclone) . the lerger-dlameter cylindrical wall 
becomes much hotter than the smaiier-dlameter one. 


<v> Turbulent combustion of ore-mixed oases within a duct 

When a pre-mlxed combustible gas burns within a duct, the 
density reduction occasioned by the exothermic reaction occasions 
a reduction in pressure with increase In distance along the duct. 
This pressure gradient causes the lighter burned gas to accelerate 
more than the heavier unburned gas. the former is therefore 
"sifted 1 * towards the duel outlet, where however the average 
temperature Is the highest. This is one of the best-documented 
examples of "counter— gradlenl diffusion": and measurements have 

clearly revealed the fact that the hot-gas fragments have higher 
velocities than the cold-gas ones (Refs 3.4. 5), 



Cvi> Sand storms 


A circumstance In which diffusion and sifting act In opposite 
directions, and maintain an approximate balance, arises when a 
strong wind blows over a sandy desert. Then turbulent diffusion of 
the Fleynofds-Prandfl kind carries sand particles from 

high-concentration regions close to the ground to 
lower- concentration regions higher up. while sifting effects a flux 
In the downward direction. 

It Is worth remarking that. In such a case, the scale of the 
fragments of upward-moving and downward— moving material Is 
much larger than that of the sand particles: the sand-taoen air 

acts like a single-phase fluid of non-uniform donslty. 

(d) A simple formulation of the slftlno orocoss 


C I) The sffttnp flux 

Let a turbulent f»u*d be taken as consisting of fragments and 
surrounding volumes of two kinds, distinguished by their having 
different densities. f>\ and P2* and different values of some scalar 
quantity ♦. 

What will now be proved Is that there Is a flux of entity ♦ in the 
direction of the local body-force vector, which Is proportional to 
the difference 4>i - 4>2* and not to the gracFenf of either quantity: 
the relative velocity U. le " ui ales enters. Here u stands for 
the absolute velocity. 

The net flux of e. viz F<$>. can be written as: 


s r)p)U)4>i * rsP2V2*2 .Cl .2-5) 

This can be re-written In terms or the relative velocity and of the 
average velocity, density and e. as follows: 


j F» = pu» ♦ <ripirgpg/p> U > | 

wherein use has been made of the definitions: 
p « r|pi + rgp2 


.O .2-6) 


.Cl .2-7) 
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p u * npyui + r 2P2 u 2 


.Cl.2-8) 


P* ■ HPl^l ♦ r 2P2&2 , (1 .2-9) 

1 * rj ♦ T 2 .Cl. 2-10> 

Equation < 1.2-6) is easy to Interpret: tne first term on the 

right-hand side represents the transfer of <t> by the average mass 
motion: and the second term represents the silting contribution. 

Further examination of that term shows that it is linearly dependent 
upon U and on (ei - • and also upon the product rirg The 

tatter, by reason of equation (1.2-10). falls to zero when either 
component is absent, as is readily understood. 

til) The relative velocity. U 

As will be explained In the next section, the magnitude of U can 
be computed by solving momentum equations *or uj and i»2 
separately, and then subtracting. However, an approximate 
formula can be obtained by neglecting a II the terms in the 
momentum equations exceot the pressure gradient, the body-force 
acceleration a. and the nterfiuld friction, as follows. The 
momentum equations reduce to: 

H ~ PI »j “ C qr 2 p llJ>U = 0 .C1.2-U1 


and ~ P 2 a] + C r\T2 P i ^ >|U = 0 .(1 .2-12) 

lex j f 


wherein C represents an inter-fluid friction factor, and 1 is a 
m e as ure of the fragment size. 


Addition of the two equations and reference to equation (1.2-7) 
lead to: 

?-Jl=0 
ex 


: ( 1 .2 13) 



7- 


wttereufion ^uintitutton In either equation l«ids to: 


U|Ui ~ C" : » <pp - p|>* t a 


.O 2-14) 


This equation Is also easy to undorstflnd: larger fragments move 

(aster, because the body force operates on the volume whereas 
the (fiction operates on the surface area; and the velocity 
difference increases with the density difference. 

the 'constant" C can be expected to Increase when the fragment 
sl/e Is small. In proportion to the reciprocal of the Reynolds 
number of the relative motion: but the assumption of constancy fs 

probably correct for hlgh-Reynolds-number flows. 

<IU> Whe n o Is line ar In ± 

Let Be stand for the sitting flux of entity #>. so that: 


Be s •fr|mm«>/ 0 ) U <42 ~ 4q ) : 0.2— 15) 

and K f>b« itoaer In e. so that: 

ip2 - Pl> ' ). .0.2-16) 

Comblrstlon of equations 0 .2-15). 0.2-16) and 
0.2-14) now yields: 

5* ' r,r ? W? - *, j 3 '* -CI.2-17* 

Since it Is also possible to relate (*2 - 4i> to 

the root -meat) square value of the fluctuations of ♦ by* 

1/2 


4' « <r|rj>> 


141 - 42» 


.0.2-16) 
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the expression for the siftlny flux in terms of 4 >' Is: 


- <r X r 2 > v '* aS |J?".j ,/a <♦* > 3'2 


(1 ? 191 


This Is a term which. In the present author's view, should niton 
appear In the deferential equation tor when this is employed In 
a turbulence model. 

Of course. It can be used only when ealimah are available oi the 
rjrg product, of the Iragmenl site t. and ot the fluctuation 
Intensity $*: but this situation Is not rare 

<lv> the sltllno Hut of mom e ntum, in Hows with 
curved stream line s 

When a fluid having velocity w follows a path having a radius ot 
curvature R. It experiences an acceleration w®/R normal to Its 
trajectory. 

If the above equations are re-worked for the situation In which the 
density Is uniform bul lire two sets of fluid fragments have d if for tug 
velocities. w| and W2. there results:^ 


Sy, * rpjifU CW2 - W|> .(!,? ?tl# 

U|*J| = C~1 (Kf2 ~ W| I I w/R .11 P PI) 

and so: 


Sw = 



fw 2 


W| 


.OP PP) 


Ot*. In terms of w* : 


S* =» (rjrg)*^ p 

tw 

1/2 

(w* 


on 

; 


iiera S w represents, of course, a 'silling shear stress'; and ii is 
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Interesting It? observe lhal Its magnitude depends upon w* while Its 
direction is always away Irom the centre of curvature, regardless 
of the direction of tlm gradient of w. 

Once again. use of the equation requires rtrg. * and w' to bo 
estimated: but Oil* Is often possible. 

2. I he mathematic al theory of the flow of two 
f otot penetrati ng contfnua 

2. 1 Background 

During recent years, it has been necessary lo predict the motion 
of steam and water In the vessels and pipes of pressurlsnd-wa»er 
systems experiencing hypothetical loss-of-coolant accidents This 
necessity has caused attention to be concentrated upon how the 
motions ol two ‘’Interpenetrating continue' can be represented 
numerically and computed efficiently: and the problem of how to 

do so has been solved. 

The llloralirre of the subject Is now too voluminous fc- summary: 
but the pioneering paper of Harlow and Amsden deserves specie 1 
mention (fief 0) : and the present author lias also made 

contributions C eg Fief 7) , 

As a consequence of these developments, ft Is possible to 
compute two-phase flow phenomena as readily as single -phase 
ones: and Indeed the ability to do so Is now embodied in a 

generalty-ayallable computer program (Fief 8.9>. which has been 
used for all the numerical work alluded lo below. 

This has been fortunate: for It has permitted the development of 

the new two-fluid turbulence model to proceed without awaiting the 
construction of a new computational toot. 

in the lollowlng section, a brief introduction to the theory will be 
provided . 

2. 2 Mathematic al f ormulation of the two-fluid model 
(a) Defin itions 

It Is supposed that two fluids. 1 and 2. share occupancy of apace 
The proportions of lime during which each can be expected lo 
occupy a particular location are the volume fractions, n and rjj. 
which obey the relation* 


r 1 * r 2 ~ 1 


. ( P 2-1) 



- 10 - 


How the (wo fluids are distinguished front one another Is . II* ary. 
however, at any location nod time, all fragments belonging to Ihtid 
1 are supposed to have the same values of tompot attire. mass 
fraction of component chemical species, velocity components, and 
fragment size. This assumption ol fluid homogeneity Is. of 
course, an Idealisation, at variance with reality: but It may well 

lie closer to the truth than do the more conventional single fluid 
models of turbulent flows. 

All the above fluid properties will be represented by the symbol <n. 
where I takes the value 1 or 9 according to the fluid. 

It will be supposed thi.t both fluids have the naum pressure at each 
point, not because ihe equation system is significantly simplified 
thereby, but because no mechanism entailing a difference ol 
pressure has been hypothesised so far. 

(bi Differential equations 

Mass balance s. There are two mass conservation laws to 
express. They can both be written In the form: 


° <p|^ q) ♦ dlV (pi 1*1 <Vi > ) = cq .4? ? 2) 

at 

wherein t stands for time. 

PI stands for the density of fluid l. 

<Vi> stands for Its velocity vector, and 

ai stands for the rate of transfer of mass Into fluid 

t from the other fluid, with which it is 

Intermingled. 

In equation (2.2-2) , instantaneous values of p. r. <V> and m are 
In question. Because *t is Impractical, in numerical calculation .. 
to compute the most rapid fluctuations ol these quaniiHo*. 
lime-average quantities are preferable. In which case lire quantity 
operated upon by dlv must be augmented by i)|. an eflncilvn 
diffusion IIuk. The equation then becomes: 


*. Cp| ri) t dlv (p( rj Vi + 0t> * sq .<? ? 3i 

at 


from which. Iiowever. over bars signifying lime averaging have boon 
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iltUMt so ms to avoid nMlie» clutter. 

« a/ at tfttm represents the dUfer«nc9 o* outflow from Inflow 
mp !«» Him dimension; dlv represents that difference Iho 

aco dimensions; and can b» rogardnd as th« negative of tho 
Inrnnro hoiwMii outflow and Inflow along an ■'Intarfluld 
nonslon* If alt these fluxes are denoted by 6(. and 0*V Is 
Hiked as a qnnoratisaUon of dlv which makes reference to aU 
air dimensions, equation (2. 2-3) can be finally expressed as: 


OIV €6| 1 - 0 .C2.2-4I 


dances of other qu antities. The above nomenclature facilitates 
o compact writing of the differential equation governing the 
sirtbMton nf the fluid property q. The equation bacemar* 


OIV (•( G| l - Ot ,f?.2-5» 


herein q h «» vniumau.c source of the entity e In phase l. and 
is the vahiA ol e associated with the this which It multiplies. 

I especial importance are the momentum equ~"~ns. In which q 
lands for a velocity component. For these equation*, the 
aim neirtc - source quantity q consists primarily of 

ies*u*o-g»adle*»t. body -for co. inter fluid- friction am) wall -friction 
arms. Thus, ter the velocity q. the equation becomes: 


□»v (u ( G 4 > * - r| «5 * bx 

ax 

* f|| Cuj — U|) — f |,w U'l .<2.2-61 

Are:- x Is the distance coordinate aligned with u|. 
p Is the prssure. 

b x Is the body- force component In the x-dl recti no. 
f(j Is the volumetric Inter — fluid-friction caef fteteot- 
U| *s ttK* local velocity of the other component . and 
f^ w Is the volumetric coefficient of friction 
bet ween fluid 1 and a nearby wait. 

•pressure 73 «adlent“ an d other kinds of diffu sion. Because of the 
association of the fluid density rj with each of the G| on the 
eft-hand side of equation (2. 2-6) . and of Its non -association with 
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the pressure-- gradient term cm to© sighs Ham! $StJ{a. -Hm effect of 
fh© pressor© gr.actferil on Pie velocity .Comp.m$m&is .of flu» iiul«to 
to inversely proportional |o stielr dwtoltlos Paossuso- gr^dter&to 
etuis© too two fluids thoielm© to have dUleatmt v^tocilteia. so tfc&a 
th©y move f©totlv©ty to 

A ,f. 

Retollv© modem intermingled fluids c ««» ^too resuft i^wn 

turbulent diffusion. to iiofn too' random iis«iv&*nant$ of MuM 

fragments. i Ho swt magnitudes of which usually propmikmai to 
the gradieni of concentration fie votum© fsaclksti) to distinguish 
UH© iv^o pftQAomma. aha? caused by to© proas***© wifi 

her© celled ^IUn9 ,t , as above* 

lc) Auxiliary rolatkma 

I rtler flui d fri ction > In oyd^ itaaf lb® ©quaUnn sHoukl fm 

~ctGb*tT. i© ?@^dared comp!®6® enough for 5uItM>on w ma9h«maiSr^ 
expressions muss h® devised for such inior- .fluid traospoti quantities 
as rrsi and Sjj. Slnc@ it to prosoot wll^r Has dtoosr^sed lives© 
qudsUoos olsawhofo Cftofs 1-0. 11, 121 - St suffices Imwo to Indicate 
what foicr.uto© bav© boon us®d to to© co^p^U&itosts ^Hicfo wiu b© 
prosocKtad botow- 

¥ho formulaikm ua®d tor fg. whlcSt appows to ©quaUtm (2 ? S> , 
to: 


U) r Cf i 1 pa IV t - Vjl .t2 2 71 

to ihh formula. eg Is i dimensionless constant. ^ is toe density 
of Ih# tighter of the two flukls: & ^ to t?m :ammmfl r >0 

Uuictof ragmen! Interlace aiea per unit voluntw of sp^rft; a«hI 
i Vt - Wjl to Uv® locai ton© average rotative spe«dd »>4 to© 
fluids 

lb© area/volum© quantity. I \ 8 s taken as toto&g pr oporumwat to 
the volume -f i action product {r$rj>. so toaS it vanishes wbon either 
tfuftl disappears. 

In more advanced work, a dfftomnUal oqsmtUm governing ft ^ to 
formulated ami solved: In I3v© present wmk However ® * has 

been taken as proparitonal to tfe® .roclptsH'al of Mae 
tllmortsion of lb© flow. 

.Snlerltuld mass transfer the Int«f-UuW mnss transfer Ins 

t>&en presumed .u proceed only to on© §Mmc»ty from 

faster ItiiiW .11 40 flower (fluid 25. This is not to© ooty possllito 
smJ mm® &ymmolflc&l owa t&v&tog 4nv*ssS6qat»d by 


-13- 


8h* pro«Mti *wn« i« oUik studios. 

rate o! »«-*■*•* ifnnstor has boon taken as directly propnrtlo'tat 
to the Inter OuM tiictkw factor. I||. too proportionality factor being 
In okc«ss of unity 

|nier!lu*d._. tvaal.. transf er and diffusion of chemical s pecies. 
Conductive beat transfer and diffusive transfer of chemical species 
between the two fluids Have boon taken as obeying a law similar to 
equation 

(d) S o lu l tao p roc edure 


The solutions which are to bo presented of the above system of 
equations have been obtained by straightforward application ©» the 
PMJFNICS computer program ( Refs S. 91 . There Is therefore no 
need Id dn more hero, by way of description of the solution 
procedure, than w* recite the main f eat ur es of what Is embodied In 
PHOEN1CS These ar«: - 

" the differential equations are replaced by fully -conservative. 

<uUy-impUcH finite domain equations, valid for a staggered grid. 

• Convection flu*es tie a‘l those represented by G| In equation 
<2.2-4T. are compuiro from ‘upwind* values of rj. pf and #t- 

• An IteraUve solution procedure Is adopted wh!ch results in 
the diminution to below pre-sot limits of all the Imbalances in alt 
the timie- domain equations. 

• A sultk lent sample of the computations fs repeated on 
successively finer grids Clo with smaller and smaller Intervals of 
space and llme> until the physically- significant results are found to 
be independent of the interval sire. 

A full account ol the computational details will be provided In 
‘PnOFNif'S I'emnnstration Reports’ from the Imperial Col'ege 
Computational Flunf Dynamics Unit. 

2. 3 Some typical results from work In progress 

(a> The plane wake 

U n Matin h»s been studying, under the author's dtrec'lon. the 
application cl the two-fluid model Of turbulence to various 
unbounded turbulent shear flows. These are dominates by mixing 
rather than siuincj: but it Is nevertheless Important to establish 

that the model wut fit the available experimental data. 

One report of this work (Ref 131 relates to the plane turbulent 
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wake, lor which experimental data are available f ip ? 3 I 
reproduced from the report, which must be referred le lor Mi 
oefaus- 

Three temperature profiles are represented In He diagram. lire 
higher represents the temperature of the hotter o* the ho fluids, 
which Is also < although (Ida Is not shown) the siowor moving cum ■ 
the tower represents tire temperature dlsiiibminr* of the 
faster-moving fluid which Is drawn Into the wake from the 
surroundings; and the curve between them represents the average 
temperature, weighted by U*e computed voiumo fractions. 

Comparison Is possible between measurements and predictions in 
respect ol the upper two curves: and It Is rather satisfactory 

The measurements relating to the upper cuivh have been marie. It 
should be mentioned, by '’conditioned sampling* 

(b) the a«J -symmetrical j et 

Fig 2.3-2. from the same report. provides cm responding 
comparisons for the avi- symmetric )et. In this case measurements 
are available for the temperature profiles of the hotter (now 
faster -moving) fluid, of the c oo l e r ( s l ower ) fluid, and of Um 
average. 

The agreement bet w e e n predictions and experiments is rather 
poorer than before; but It Is premature to attempt to ascribe a 
cause. 


to) A combustion study 

"I he two-fluid model of turbulence Is especially welt suited to the 
representation of combiisllr ■ processes: tor It permits explicit 

allowance to be made feu the fact that damns consist of 
Intersperslons of holler and colder fragments til gas and that thr* 
chemical reaction rate Is especially active at tire interfaces between 
fragments 

Figs 2. 3-4. 5. 6. 7. 8 and O contain some rostilis ol a computation 
ol whal may happen when a shock wave.’ Induced by the 
Obstructing of flow In a pipe, travels through an Inlersporsfon of 
hot-gas fragments In a surrounding colder combustible gas fig 
2.3-3 sketches the apparatus In question. 

Although reference must be made to a more complete report t-ftet 
14) for a full explanation. It can perhaps be perceived, by 
Inspection of the i J, Bgrams and their legends, that: 

• a pr* 'sure wave ol growing ampHturfe paves along tie 
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ducl: 

• thfl volorlty difference between the hot end cow qas 
inciaaw resulting In e 'sliltog* of hot gee Into the 
power - pcA^uia region: 

e *t is in the region where the velocity difference is a r notes t 
thnt the inter fluid mess transfer, and the reaction rate, are 
most vigorous 

CO) Other studies 

The a*‘itior‘s student* and. co-workers are pursuing further studies 
in this area. Attention Is being paid to: 

e turbulent Hows near walls: 

e turbulent dlffitsir. » flames under tree-convection conditions: 

e IV»ws In hydrocyefones: 

• Users Induced by Rayleigh— Taylor Instabilities. 

fV'porls so this work are now beginning to appear. In lite first 
Instance In the Imperial College CFDU Report series. 

f 

Ti » work proceeds slowly, because the aim Is to develop » single 
tw -fluid model which is capable of titling a large number of 
dh use experimental data. However, progress h being made. 

3. LOME JJPEN. QU r j.IONS 

The tw-*-l» rid model of turbulent (lows Is not yet ready for 
recommendation to practising eng'neers as a reliable predictive 
too': fr>«. although plausible-seeming predictions have been 

made. and qualitative agreement with experiment has been 

achieved in respect r»t features which other models cannot predict 
at all •*: ere are many questions still to be answered. Some of 
the qves.tons to which the present author attaches Importance will 
here be flsied. They are: — 

Ca) What is the mo-t s titablv formulation for the Inter-llutd 
mars 'ransfer taw l seems certain that early presumptions 
of on* i way transit*. C eg from slower to fester, nr from 
low-ve-*Uclty * ■ hlgh-vortlclfy) cannot be universally 
sattsiacic— ut should the mass-transfer rale c nance sign 
when the •luma fraction of a fluid attains 0. 5. or some 
other number dependent upon density ratio or other 
dimensionless quantity? 

<b> "'Hi It suffice to employ a single length scale which Is 
characteristic of both fluids? And. If so. ht>w shmitd the 
lengt t icale be formulated to allow for the experimental 
observations that larps fragments 'swallow* small ones when 
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Ihs silting process Is dominant* and also In some shear r 
(lows («9 tl»o piano mining layer* . whereas at vary high 
QbMi ratoa ths average fragment sirs dlnitnlsties 

Although it seams wise to develop the model try use ol the 
full set ol two- fluid equations, as dasr-ilhed hr section 2 9 . 

It would bo very desirable to provide a model employing only 
pie conventional single lluM equations with additional tar ms 
tor the sifting tlur II this is to be dona, whai means 
should be used tor calculating etc. which appear In the 
expressions for this flux? 

td* INhal experimental data already exist* ami what new 

experiments can usefully be devised, which will enable lire 
as- yet* unknown constants and functions of the two field 
model to be most easily and accurately deduced? 

(e* What rote can computational fluid dynamics (day. when 

applied to fragment formation, -distortion and Inter sr linn 

processes, to elucidate the unknowns Irr the iwu fluid model 
equations? And how can such a program ol study best be 
formulated - and funded? 

(II Will it indeed be possible to advance to the point of 

providing engineers with a valuable design t«mt without 
Increasing the number ol fluids considered beyond two? 
And. If not. what Is the smallest numlter ol fluids which will 
euillce? 

4. CONCtUOtNQ REMARKS 

It Is the present author's opinion thal. despite the slowness of the 
progress and the number of still unanswered questions. Ua 
tWQ-Uuld Idea can form a good basis for further advances 

As to slowness of progress. It must be observed that 
turbulence-model development at the single liutd level has 
stagnated In recent years. An increasing number of -experimental 
studies is showing the loedequacy Ol the k «* mode! as formulated 
In the early 1970 s: yet generally valid refinements ami 

Improvements have mil been forthcoming; and thn 

Reynolds- stress* modelling approach. In tho writer's opinion. 
Increases the difficulties of making predictions wtiimut significantly 
Improving their qualify 

Because ahllily lo perform two- fluid raiculallmis has. until 
recently, been con lined to a tew persons and Institutions heavily 
engaged In other work, the nan-poww devoted to pursuing 
Iwo-tluid-modeUUig has remained small. rktw tmwever u Is worth 
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nranHonlng. Ih« Iron availability of PHOENICS to academic 
Institutions makes II possible for many more individuals and loams 
to Join In the tnplctttlon of this now branch of turbulence 
rosoarch. should they wish to do so. 
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Flgure 2. 3-3. 

(a) Sketch of a tube containing an Intersperslon of hoi 

combustion products In cooler pre-mlxed combustible gas. 
Initially both gases flow at equal velocity from felt to right; 
then a valve at the right-hand end closes, senoing a 
pressure wave along the duct. Combustion is thereby 
Intensified. 

<b> Sketch of the pressure and combustion wave trajectories on 
a distance-time plane. 

<c> Sketch of Instantaneous profiles of pressure p. Interfluid 

mass-transfer rate m. cold- gas vetcity uj. and hot-gas 
velocity U2* 1 ms alter dosing of the valve. 
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<F*guro 2.3-4, The flame of Fig 2.3-3* 0. 2 milliseconds after the 
closure of the right-hand end. 
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Figure 2.3-5. The flame of Fig 2.3-3* 0.4 milliseconds after 
closure of the right-hand end. 
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Figure 2.3-6. The flame of Fig 2. 3-3. 0.6 milliseconds after 
closure of the right-hand end. 



Figure 2.3-7 The (lame of Fig 2.3-3. 0 8 milliseconds after 
closure of the right-hand end. 
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Figure 2.3-8. Distributions of pressure, velocities . the 

mass-transfer rat® and other quantities, 1 millisecond alter 
right-hand-end closure. 
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Figure 2.3-9. Variations with time o' pressure, reactedness. 
mass-transfer rate, volume fraction ol hotter fluid (V) , velocities 
of hotter <H) and colder <Q fluids, all for a location 3/4 ol the 
distance along the pipe, ie 25 cm front the closed end. 



